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Abstract. We study the effect of time-odd components of the Skyrme energy density 
functionals on the ground state of finite nuclei and in nuclear matter. The spin-density 
dependent terms, which have been recently proposed as an extension of the standard 
Skyrme interaction, are shown to change the total binding energy of odd-nuclei by only 
few tenths of keV, while the time-odd components of standard Skyrme interactions give 
an effect that is larger by one order of magnitude. The HFB-17 mass formula based 
on a Skyrme parametrization is adjusted including the new spin-density dependent 
terms. A comprehensive study of binding energies in the whole mass table of 2149 
nuclei gives a root mean square (rms) deviation of 0.575 MeV between experimental 
data and the calculated results, which is slightly better than the original HFB-17 mass 
formula. From the analysis of the spin instabilities of nuclear matter, restrictions on the 
parameters governing the spin-density dependent terms are evaluated. We conclude 
that with the extended Skyrme interaction, the Landau parameters Gq and G'q could 
be tuned with a large flexibility without changing the ground-state properties in nuclei 
and in nuclear matter. 
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1. Introduction 

Despite many theoretical and experimental investigations, the spin and the spin-isospin 
channels in either the ground and the excited states of nuclei are still widely open for 
future study [U [21 El IH El [6] . It is indeed difficult to probe the spin and the spin-isospin 
channels of nuclear interaction since the ground states of nuclei are non-spin polarized 
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in the case of even-even nuclei and at most polarized by the last unpaired nucleons in 
odd nuclei. 

The analysis of spin and spin-isospin collective modes such as magnetic dipole 
(Ml) and Gamow- Teller (GT) states gives access to the nuclear interaction in these 
channels. The Landau parameter Gq has been deduced from the analysis of the GT 
mode: a model based on Woods-Saxon single-particle states plus one-pion and rho 
meson exchange interactions gives Gq = 1.3 ±0.2 (see Ref. [21 [3] and references therein). 
A slightly different value Gq = 1.0 ±0.1 was derived from observed GT and Ml strength 
distributions using the phenomenological energy density functionals DF3 [1, 7J. Anyway, 
in both cases empirical single-particle energies are used (that is, for states close to the 
Fermi energy the effective mass m*/m is ~ 1). 

Self-consistent Hartree-Fock (HF) plus Random Phase Approximation (RPA) 
calculations constitute a somewhat different framework, in which the density of states 
around the Fermi energy is lower (or, equivalently, the effective mass is about 30% 
smaller). As compared with the empirical case, the unperturbed particle- hole transitions 
have larger energies and as a consequence one needs a smaller residual repulsive effect to 
fit the observed GT peak. It is not surprising, therefore, that self-consistent calculations 
of the GT resonance, performed using different Skyrme interactions in Ref. [6], point 
to Gq ~ 0.6. Skyrme interactions are characterized by a spin-isospin G[ parameter as 
well, and specific terms of the effective mean field, like the spin-orbit potential, may 
also break the simple correlation between the GT properties and the parameter Gq. 
However, as widely used interaction like SLy5 have unrealistic (negative) values of Gq, 
it is undeniable that adding more flexibility to the spin-isospin part of Skyrme forces is 
useful. 

The spin and spin-isospin component of the nuclear interaction is also reflected 
into the time-odd component of the mean field. The properties of even nuclei give 
constraints to the time-even component of the mean field while very little is known 
about properties of the time-odd mean fields. Time-odd components of the mean field 
compete with pairing correlations in determining the odd-even mass staggering [8]. Fast 
rotation induces time-odd components in the mean field [9] which could be probed from 
the measurement of the dynamical moments of superdeformed bands. 

In Ref. [lOj, new spin-density dependent terms have been introduced on top of 
standard Skyrme forces in order to remove the ferromagnetic instability associated with 
all Skyrme parameterizations [5J. The new terms retain the simplicity and the good 
properties of Skyrme interactions for nuclear matter and the ground states of even- 
even nuclei. However, these new terms slightly change the properties of odd systems. 
This work aims at studying quantitatively the effects of these new terms on the ground 
state properties of odd nuclei, in particular the total binding energy and the density 
distribution. Since these terms contribute only in odd nuclei, both odd-even and odd- 
odd nuclei are considered in the present study. To provide an approximate maximal 
estimate of the effects while keeping our model simple, we perform HF calculations with 
the following approximation to treat odd nuclei. We use the equal filling approximation. 
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so that the time-reversal symmetry is not broken, but with an additional ansatz. 
As within this scheme the spin-densities would be by definition equal to zero, when 
constructing the spin-densities with the wave function of the odd nucleon, we assume 
that the spin-up state is completely filled while the spin-down state is empty between 
the two possible spin orientations (or, equivalently, the opposite). We call this procedure 
for the construction of spin-densities the one-spin polarized approximation (OSPA). The 
OSPA gives an upper value of the contribution of the new terms. 

The article is organized as follows: in Sec. II we remind the newly proposed spin- 
density dependent terms as in Ref. [10]. In Sec. Ill, we estimate how much these new 
spin-density dependent terms affect the total binding energy of odd nuclei. In Sec. IV, 
it will be shown that, introducing these terms in the most predictive HF-Bogoliubov 
(HFB) mass formula |Tl], the quality of the mass fit can be recovered with an optimal 
renormalization of the Skyrme parameters. In Sec. V we will analyze the ground state 
properties of infinite nuclear matter with respect to spin-polarization and give a range 
for the parameter of newly introduced spin-density Xg. Finally, conclusions and outlook 
are given in Sec. VI. 

2. Extended Skyrme interactions 

As described in Ref. [lOj, the new spin-density dependent terms added to the 
conventional Skyrme force are of the following form: 



where Pq- = {l + ai-a2)/2 is the spin-exchange operator, r = ri — r2 and R = (ri + r2)/2. 
In Eq. ([T]), we have introduced the spin-density ps = Pt ~Pi spin-isospin-density 
Pst = Pn] — Pni — Pp] + Ppi- Spin symmetry is satisfied if the power of the density- 
dependent terms 7^ and 7^^ is even. 

The total energy -Etot in finite nuclei is related to the (local) energy density 'H{r) 
through 



In the following, we adopt the notation of Ref. [12] where Ti is expressed as the sum 
of a kinetic term /C, a zero-range term Tig, a density-dependent term 7-^3, an effective- 
mass term TieS, a finite-range term T^sn, spin-orbit and spin- gradient terms (7-^so and 
TYsg), and eventually the Coulomb term TYcoui- However, the expression for the energy 
density provided in [12j holds in the case of time-reversal symmetry; in odd nuclei, the 
energy density Ti acquires also a dependence on the spin-densities ps and pst, hereafter 
named 7i°'^*^, even without additional terms in the force depending on these densities. 
Expressions for 7i°'^'^ have been derived for instance in Refs. [9l [I3] or in Appendix I of 
Ref. [Hj. For the reader's convenience we repeat here, in Appendix A, the expression 



y^'^^ri.r^) = \ti{l + xlP,)[pmr5{v) + ^^^(1 + 4*^<.)[P.t(R)]^"^(r) 


(1) 




(2) 
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The additional terms ([T]) modify the standard H3 contribution to be 0.3 + 
where the last two terms read 
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where q = n,p. 

As already mentioned, the additional contributions to the mean field are zero in 
even-even nuclei. Since most of the Skyrme interactions are adjusted on (few) even- 
even nuclei, it is thus possible to add for these interactions the new terms ([1]) in a 



perturbative manner. The new four parameters t 



3; ^3; ''3 



ti^ and X3* in Eq. ([T]) have been 



adjusted in Ref. [T^ in order to reproduce the Landau parameters extracted from a G- 
matrix calculation in uniform matter, while 7^ = = 2 is imposed by spin symmetry. 

Some other interactions, in particular those produced by the Brussels-Montreal 
group, are globally adjusted by fitting the properties (essentially the nuclear masses) of 
both even and odd nuclei. In this case, when including the new spin-density terms, a 
global re-adjustment of the interaction might be necessary. 

In the following, we first analyze the effect of the new spin-density dependent 
terms for a few selected nuclei using the SLy5 Skyrme interaction for which the new 
spin-density dependent terms are added perturbatively. Later, on the basis of the 
latest BSkl7 Skyrme parameter set, we study the global impact of our new terms on 
nuclear masses and show that the Skyrme parameters can be refitted to provide almost 
equivalent properties. 



3. Ground states of Ca nuclei 



According to the OSPA, we define the densities Ps(r) and Pst{r) in odd nuclei as 

Psir) = , (6) 

i 

Pst{r) = Yl ¥^Kr)ms{i)mt{i), (7) 

i 

where ms{i) and mt{i) are the spin and isospin z-component for each single nucleon 
having the wave function (pi{r). The last occupied state fully contributes to the spin- 
density. It is then clear that the OSPA corresponds to maximizing the spin-density and 
its effects. 

Note that if the time-reversal symmetry is not broken (in the filling approximation, 
for instance) both spin-up and spin-down states must be degenerate and the densities 
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Ps and Pst are zero. Since in this work we aim at an approximate and maximal estimate 
and not at a precise prediction of the effects, we use the HF method with the OSPA to 
treat odd nuclei. 

All calculations in this Section are performed with the SLy5 parameterization for 
the conventional Skyrme force and the additional parameters for the spin-density terms 
given in Table 2 of |10j. We choose two systems: ^^Ca and ^^Ca. Being built on top 
of double magic nuclei, pairing is neglected for these nuclei in this section. Both nuclei 
are odd-even, then the proton spin-density is zero and the spin-isospin-density is equal 
to the spin-density for neutrons. For the nucleus ^^Ca the spin-density is built with 
the neutron I/7/2 wave function while in "^^Ca it is constructed with the neutron 2^3/2 
wave function. The choice of these nuclei has been driven by the fact that one of the 
nuclei, ^^Ca, has a spin-density which probes the external region of the nucleus, while 
for the other nucleus, ^^Ca, the spin-density may probe more the central part. We could 
not use the OSPA if the single-particle wave functions changed when adding one or two 
nucleons. Therefore, as a necessary step to proceed with the evaluation of the impact 
of spin-density dependent terms within the OSPA, we have first checked that the wave 
function of the neutron state I/7/2 does not change appreciably when passing from the 
even nucleus ^''Ca to the next even isotope with + 2 neutrons, ^^Ca. An analogous 
check has been made for the wave function of the proton state I/7/2 in ^°Ca and in 
the Z + 2 isotone ^^Ti. Also, the spin-density calculated with the OSPA for "^^Ca has 
been compared with half the difference of the neutron densities of ^^Ca and ^°Ca, where 
effects coming from the rearrangement of the deeper state are also included. Very small 
differences are found, mainly in the central region. Analogous results are obtained for 
^^Ca, ^^Ca and ^^Ca. 

The spin-densities in ^^Ca and ^^Ca are plotted in the two panels of Fig. [1] together 
with the neutron densities of both nuclei. For comparison, the neutron densities of 
the nearest even-even nuclei ^^Ca and ^°Ca are also displayed. The spin-densities 
and the neutron densities have been calculated in the two odd-even systems. Two 
independent calculations have been actually performed with and without the new spin- 
density dependent terms ([1]) and only negligible differences have been found so that they 
are not appreciable in the figures. 

Since the contributions to the total energy Etot coming from the spin-density 
dependent terms ([3]) and (jlj) are proportional to the square of the spin densities 
represented in Fig. [H these contributions are expected to be negligible as compared 
to the usual density dependent term H-^. To be more quantitative, we have made 
several calculations for the ground-state energies that are summarized in Tables [T] 
and [2j The partial contributions to the ground-state energy (j2]) are written Emf = 
J d^r {Ho + ^3 + n,s + ^fin + ^sg), £^so = J d^r H^o, ^coui = J d^r Hcoui and 
E\^in = J d^r Tikin, using the notations of Ref. [T2]. In Table [H the total energy, the mean 
field, the spin-orbit, the Coulomb and the kinetic contributions to the total energy and 
the single-particle energy (for the neutron states I/7/2 or 2^3/2) are provided for the odd 
nuclei ^^Ca and '^^Ca and for the nearest even-even nuclei *^^Ca and ^°Ca. The calculation 
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for the odd nuclei are performed either with or without the spin-density dependent 
interaction ([1]). The difference between the total energy, and its contributions, without 
and with the corrections AE is shown to be less than 50 keV in both odd nuclei ^^Ca 
and ^^Ca. It is interesting to notice that the kinetic energy contributes to reduce the 
impact of the spin-dependent interaction ([T]). The spin-orbit and the Coulomb energies 
are very weakly affected by the spin-density terms. 

In Table [2] the total and separate contributions coming from the time-odd terms 
of the Skyrme interaction (see Eqs. ( ]A.2[) -( lAl5l) in Appendix A) are provided, in the 
case of the odd nuclei ^^Ca and "^^Ca: these contributions are calculated perturbatively 
within the OSPA. These terms are classified according to the standard notations and 
labeled as Eq'^'^, -Ecg+fin+sg ^3'^'^ which are the contributions from the central, the 
momentum dependent and the density dependent terms, respectively. The corrections 
Eq'^'^ and E^'^'^ give a dominant and repulsive contribution which increases the total 
energy while the correction -Eeff+fin+sg smaller and attractive. The total correction 
remains quite small, that is, of the order of 0.15-0.3 MeV for both nuclei. Note that 
the sign of these corrections could change from one Skyrme interaction to another, but 
such corrections in Table [2] remain larger than Ai^TOT in Table [H From the quantitative 
comparison shown in Tables [1] and [2], we can infer, as expected, that the new spin-density 
dependent terms ([1]) modify the ground state energies of odd nuclei much smaller than 
those coming from the time-odd terms of the standard Skyrme interaction. 



4. Global adjustment on the nuclear chart 

Some Skyrme interactions have been determined by fitting the parameters to essentially 
all of the available mass data and therefore are constrained to even as well as odd 
systems. In this case, the new spin-density dependent terms added to the standard 
Skyrme interaction may modify the quality of the fit. To study the impact of the new 
terms on the prediction of nuclear masses, we consider now the latest and most accurate 
HFB-17 mass formula (with a rms deviation of 0.581 MeV on the 2149 measured masses 
of [TTj) obtained with the BSkl7 Skyrme force |llj . 

If we consider the parameters of the additional spin-density dependent 
interaction determined in Ref. [lO], namely = 2 x 10^ MeV fm"^, tf = 1.5 x 
10^ MeV fm^, X3 = — 2, X3* = and 7^ = = 2, we find that the impact of the new 
terms on nuclear masses are relatively small, as already discussed in Sect. III. Fig. [2] 
shows the mass difference obtained by a spherical HFB calculation when the spin terms 
are added or not. The nuclear masses of odd- A and odd-odd nuclei are globally increased 
by a value of the order of 100 keV. For light nuclei this correction is the largest and 
can reach at most 350 keV. As already shown in Ref. |inj, the spin-density terms are 
repulsive and leads to an increase of the rms deviation with respect to all the 2149 
measured masses from 0.581 MeV to 0.591 MeV, keeping the good quality of the mass 
fit. Deterioration can potentially be avoided if the force parameters are re-adjusted to 
the nuclear mass data. 
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We have determined a new BSkl7st Skyrme force which essentially corresponds to 
the BSkl7 force, but for which the Skyrme as well as the pairing parameters have been 
slightly renormalized by a new fit on the whole set of mass data. The new parameter 
set has been built with the idea of evaluating a maximum effect of the new spin-density 
dependent term (in the same spirit of the OSPA discussed above). The values adopted 
for the parameters are ^3 = 4 x 10"^ MeV fm^, tg* = 3 x lO'' MeV fm^, = —0.5, 
xf = and 7^ = 'jst = 2. ^3 and tf are twice larger than in Ref. [10] and this leads 
to Landau parameters Go = —0.03 and G'q = 0.99 in symmetric matter at saturation 
density, whereas the parameter set of Ref. [10] is associated with lower values, namely 
Go = —0.36 and Gq = 0.75. As explained in the Introduction, the value of Gq may 
be considered large in keeping with the effective mass 0.8 of BSkl7, yet still quite 
acceptable for the purpose of the present study. With a value of twice larger than 
BSklGst, the value of Xg = —2 in Ref. [10] has been modified consistently to be —0.5 
to keep the contribution to the Landau parameter Gq in spin-saturated infinite neutron 
matter identical to the one determined in Ref. [TU] [see Eq. (11) of Ref. [ID]]- As far 
as the parameter 0:3* is concerned, there is so far no constraint that could guide us in 
fixing its value. For this reason the zero value was assumed in Ref. [TU]. However, as 
shown in Sec V, the value of xf influences the stability of the partially polarized neutron 
matter, i.e., the lower its value, the higher the barrier between the S = and S = 1 
configuration. For this reason, as discussed in Sec. V, the value of is set to —3. 

The strategy of the mass fit is the same as the one described in Ref. [TT| [T6] . 
In particular, the mass model is given from deformed HFB and the pairing force 
is constructed from the microscopic pairing gaps of symmetric nuclear matter and 
neutron matter calculated from realistic two- and three-body forces, with medium- 
polarization effects included. To accelerate the fit, a first estimation of the energy 
gained by deformation is performed and subtracted to the experimental masses. A 
first series of parameters are then obtained from the comparison of spherical HFB mass 
model with the corrected experimental masses. The corrections due to deformation 
are then reevaluated and a new fit is performed. This fast procedure is repeated until 
convergence. The isoscalar effective mass m*/m is constrained to 0.80 and the symmetry 
energy at saturation J to be 30 MeV in order to reproduce at best the energy-density 
curve of neutron matter [TJ] from realistic two- and three-nucleon forces. Note that the 
parameters of the additional spin-density interaction are not fitted in this procedure. 

The final force parameters labeled BSkl7st and resulting from a fit to essentially all 
mass data are given in Table [31 It can be seen that there is little difference between the 
parameters of the BSkl7 and BSkl7st forces; note that the parameters of the rotational 
and vibrational corrections [16], are identical for both forces. The same holds for the 
parameters of infinite matter with the incompressibility coefficient = 241.7 MeV, 
the volume energy coefficient = —16.053 MeV and the isovector effective mass of 
m*Jm = 0.784, as in Ref. [II]. 

The rms residuals for the BSkl7 and BSkl7st sets are compared in Table HI The 
inclusion of the new spin-density dependent terms which slightly deteriorated the 
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accuracy of the BSkl7 force leads now even to a small improvement of the predictions 
by about 6 keV, with the BSkl7st force parameters. Both forces have been used to 
estimate the mass of the 8508 nuclei with 8 < Z < 110 and lying between the proton 
and neutron drip lines. Differences of no more than roughly ±0.5 MeV are found on the 
entire set. 

Finally note that the nuclear binding energies remain extremely insensitive to the 
value adopted for the parameters and xf. In particular, setting = —2 instead of 
—0.5 decreases the binding energy by no more than 10 keV. Similarly, a change of xf 
from zero to —3 impact the masses by maximum 8 keV. Therefore, these two parameters 
should be constrained rather by stability conditions of polarized or non-polarized infinite 
nuclear matter, as shown in the next Section. 

5. Ground state of infinite nuclear matter 

In Ref. [To], it has been shown that the new spin-density dependent interaction ([1]) 
stabilize non-polarized matter with respect to spin-fluctuations. As shown in Fig. [3l 
with the new terms, the Landau parameters Gq and Gq remains larger than —1 at 
all densities for SLy5st, LNSst [19] (which include the new terms as parametrized in 
Ref. [TO]) and the BSkl7st forces 

However, it has not been checked if the true ground state is really that of non- 
polarized matter. To do so, the energy for different spin-polarizations should be 
compared with that of non-polarized matter. This is done in the next two subsections 
for both symmetric nuclear matter and neutron matter 

5.1. Symmetric nuclear matter 

The difference of the binding energy of spin-polarized matter to that of spin-symmetric 
matter, E/A{6s,p) — E/A{6s = 0,p), is represented in Fig. H] as a function of the 
polarization 6s = (pt ^ Pi)/ P- the left panel, we have represented the binding energy 
of the BSkl7 Skyrme interaction without the new spin-density dependent terms ([T]). 
The instability occurs between p=0.18 and 0.2 fm~'^. At p=0.2 fm~'^ the minimum 
energy is obtained for a polarization 55=0.76. 

The binding energy of BSkl7st which includes the spin-density dependent terms 
is represented in the right panel of Fig. HI As expected, the energy of non-polarized 
matter is convex around ^5 = 0, but there is a change of convexity for large values 
of 6s ~ 0.8. We have indeed observed a large influence of the parameter Xg on the 
binding energy of fully polarized matter. In Fig. [5] are represented the binding energies 
E /A{5s, p) — E /A{5s = 0, p) for the three modified Skyrme interactions BSkl7st, LNSst 
and SLy5st for which we changed the values of the parameter x%,. Its values are indicated 
in the legend of Fig. [5l We fixed the density p=0.6 fm~^ to be the highest value where the 
nuclear Skyrme interaction is applied. For values of the parameter a;3=-3 the ground 
state of nuclear matter is fully polarized {5s = 1) for the interactions BSkl7st and 
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LNSst. Increasing the value of the parameter Xg from -3 to 0, the binding energy of 
fully polarized matter is going up in Fig. [51 There is then a critical value above which 
non-polarized matter is the ground-state of nuclear matter. 

One could obtain an estimate of this critical value by analyzing the contribution of 
the new spin-density terms ([T]) in spin-polarized symmetric matter. It reads 



t 



2xi - 1 



Kisym.) = ^p'pf [1 + ^^Si] , (8) 

and 1-L'^{sym.) = since pst = 0. The term ([8]) is zero for the spin-symmetric matter 
with ps = and is always positive for = 1 if one chooses Xg > — 1. It is thus clear 
that one necessary condition for the spin-symmetric matter to be the absolute ground 
state at all densities is Xg > —1. This is the case for the adopted value of Xg = —0.5 
for BSkl7st. Nevertheless, as shown in Fig. [5] for instance for SLy5st, the stability of 
spin-symmetric matter could be obtained even if Xg < —1 at the density p=0.6 fm~'^. 
At lower density, SLy5st is also stable, but not at higher density. 

We remind that from the analysis of the Landau parameters the stability around 
spin-symmetric matter requires that Xg < 1 (see Eq. (11) of Ref. [lO]). As a conclusion, 
one could adjust the parameter Xg inside the range — 1 < Xg < 1. 

5.2. Neutron matter 

The case of pure neutron matter is somehow very peculiar. The correction due to the 
spin-density dependent terms reads 

ni{neut.)=^lp'pf{l-xl)[l-6l\. (9) 

It is then clear that the correction is zero for = and also for ^5 = 1. This property is 
related to the anti-symmetrization of the interacting nucleons. Indeed, in fully polarized 
neutron matter, the quantum numbers for spin and isospin are S = 1 and T = 1 while 
the new spin-density dependent interaction ([T]) act in the L = channel. The new 
spin-density dependent interaction ([1]) have thus no effect at all in the purely spin- 
polarized neutron matter. Only odd L terms could play a role in the fully polarized 
neutron matter. This property has been used to provide a necessary condition to remove 
the spin instabilities and lead to the condition —5/4 < X2 < —1 [20]. This condition 
has been used in the fitting procedure of SLy5 [12] and it explains the robustness of 
the spin-symmetric ground state for this interaction. However if more flexibility in the 
Skyrme parameters is necessary, it might be interesting to introduce an interaction of 
the following form 

t^(l + x^P.)k'p.(R)-5(r)k. (10) 

This L = 1 term will not contribute to spin-symmetric matter and could be adjusted to 
fit the energy of fully polarized matter. 

In Fig. E], it is shown that the new spin-density terms ([T]) contribute to the binding 
energy for partially polarized matter and tend to stabilize the state 6s = 0. After the 
ferromagnetic transition, the state ^5 = is not any more the absolute ground state in 
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pure neutron matter. However, the new spin-density terms generate a potential barrier 
between the non-polarized and fully polarized states. The height of the barrier depends 
on the chosen parameters of the new spin terms ([T]) as well as on the neutron matter 
density. In the right panel of Fig. [6l the top of the barrier at the ferromagnetic transition 
is located around the polarization 6s = l/\/2 with a height (in MeV per nucleon) of 

^{tl{l-xl)+f,\l-xf))p' . (11) 

So, the higher the density, the higher the barrier. This barrier height is always positive 
since the curvature of the binding energy around ^5 = is related to the Landau 
parameter Gq^nm in neutron matter which is larger than -1 at all densities (see Fig. [3]). 
For the BSkl7st force, at p ~ = 0.19 fm~^ this barrier amounts to about 10 MeV 
per nucleon. 

From Eq. (ITTi) . it can also be seen that the parameter ^3* influences the height of 
the barrier. The lower xf, the larger the barrier. For BSkl7st, we set xf = —3 to 
get a barrier above the non-polarized ground state of the order of 10 MeV per nucleon. 
This condition is chosen with respect to the theoretical predictions [2n [22l [23| [23] 
that nuclear matter is spin-symmetric up to reasonable high densities (see for instance 
discussion in the introduction of Ref. ^Q\)- With a barrier height of the order of 10 MeV 
per nucleon, newly born neutron stars with typical temperature going from 1 to 5 MeV 
might not be the site of a ferromagnetic phase transition. The transition towards a 
non-polarized cold neutron star shall then be stable and the remaining neutron star 
spin-symmetric. Notice however that despite the theoretical predictions that dense 
matter is not spin-polarized [211 ISSl [23l [23] , there are no strong evidences against the 
occurrence of ferromagnetic phase transition from observation of neutron stars. Indeed, 
a spin-polarized phase in the core of neutron stars might induce the very huge magnetic 
fields 10^^^^^ G yet unexplained that have been proposed as the driving force for the 
braking of magnet ars [25] . 

6. Conclusions 

The occurrence of the spin instability beyond the saturation density is a common feature 
shared by different effective mean-field approaches such as Skyrme HF, Gogny HF [26] 
or relativistic HF [27] . The analysis of the spin component of the Skyrme interaction as 
well as its extensions might thus guide us to a wider understanding of the spin channel in 
general for nuclear interaction. There are many reason for looking at this channel. For 
instance for its competition with pairing correlations in the odd-even- mass staggering [8], 
for rotating superdeformed nuclei [9], for a better description of GT response, and for all 
applications in astrophysics such as for instance predictions of /3-decay half-lives of very 
neutron rich nuclei produced during the r-process nucleosynthesis [7], rehable calculation 
of neutron star crust properties such as ground-states and collective motion [28], for Ou— 
and 2z/ double beta decay processes, for URCA fast cooling, and also for neutrino mean 
free path in proto-neutron stars. 
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In this paper we have carefully analyzed the ground state properties of finite nuclei 
and infinite matter obtained by the extended Skyrme interactions with the spin-density 
terms proposed in Ref. [10]. In finite even nuclei, the new spin-density interaction ([T]) 
are simply zero and from the OSPA we have shown that these terms has only negligible 
contributions to the ground-state of odd nuclei. These results has been obtained either 
by introducing the new spin-density dependent terms in a perturbative way to existing 
Skyrme interactions such as SLySst or performing a global adjustment of the parameter 
set on the nuclear chart. A new mass formula HFB-17st adjusted in the whole isotope 
chart (2149 nuclei) is obtained with the rms deviation of about 575 keV. From the 
analysis of the ground state of nuclear matter, a range for the parameter x\ is restricted 
to — 1 < X3 < 1 in order to stabilize the spin-symmetric matter. The case of neutron 
matter is also discussed and it is shown that the new terms ([T]) with the relative angular 
momentum L = have no contribution to fully polarized neutron matter. Thus, it has 
been shown that by using the extended Skyrme interactions [10] , the Landau parameters 
Go and Gq could be tuned to realistic values without altering the ground-state properties 
in odd nuclei as well as of nuclear matter. 
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Appendix A. Time-odd components in the mean field of the Skyrme 
interaction 

In odd nuclei, the energy density 7i acquires a dependence on the spin-densities ps and 
pst [HI [13]. Respecting the decomposition of the Skyrme energy functional proposed in 
Ref. [12], the components of Ti"'^^ are 



-T/odd 



1 

4«° 



■odd 



u 



-T/odd 



24 ^'^ 



(A.l) 
(A.2) 



-T/odd 



-^[-tl(l - 2Xi) + ^2(1 + '2X2)]PsTs + ^[^2 - ti]p,tTst , 



l-tl(l - 2Xi) + t2(l + '2X2)\PsTs + 4, 

lb lb 

= - 2x,) - t2(i + 2xM + - t2]j: 

where (r^t) is the spin (spin-isospin) kinetic density energy defined as = — 
{Tst = T„| — — Tp^ + Tpj) and (j^J is the spin (spin-isospin) current. 



\2 



(A.3) 
(A.4) 
(A.5) 
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Table 1. Total energy, mean field, spin-orbit. Coulomb and kinetic contributions to 
the total energy (third column) and single-particle energy of the neutron state I/7/2 
for "^^^^^Ca and 2^3/2 for "'^"^"Ca calculated, for the nearest even nuclei and for the 
odd nuclei, without /with the spin-dependent terms ([1]) in the mean field. AE is the 
difference of energy with and without the spin-dependent terms (HI). 



Nucleus 


Etot 


Emf 


Eso 




-Ekin 


s.p. energy 




(MeV) 


(MeV) 


(MeV) 


(MeV) 


(MeV) 


(MeV) 


42 Ca 


-362.591 


-1111.434 


-9.173 


72.023 


685.993 


-9.66 


4iCa 


-352.942 


-1081.395 


-5.259 


72.116 


661.596 


-9.64 


41 Ca with (P) 


-352.918 


-1081.359 


-5.259 


72.115 


661.584 


-9.64 


AE 


0.024 


0.036 


0.000 


-0.001 


-0.012 




50Ca 


-429.654 


-1326.381 


-33.958 


70.905 


859.779 


-5.84 


49Ca 


-423.876 


-1305.865 


-33.639 


71.105 


844.523 


-5.70 


49 Ca with (P) 


-423.825 


-1305.754 


-33.634 


71.102 


844.461 


-5.70 


AE 


0.051 


0.111 


0.005 


-0.003 


-0.062 





Table 2. Total and separate contributions to the energy from the time-odd (spin 
symmetry breaking) terms of the SLy5 Skyrme interaction |14l 113] . 



Nucleus 


Tpodd 
-^TOT 

(MeV) 


(MeV) 


Tpodd 

-^off-|-fin4-sg 

(MeV) 


^odd 

(MeV) 


44 Ca 

49Ca 


0.329 
0.151 


0.196 
0.187 


-0.007 
-0.176 


0.140 
0.140 
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Table 3. Parameter sets for BSkl7 and BSklTst: the first 15 lines give the Skyrme 
and additional spin parameters, the 16th the spin-orbit term, the 17th to 20th lines 
the pairing parameters and the last 4 lines the Wigner correction [16) . See the text 
and Ref. [11 for more details 







JjoKl ( Sl 


to 


-1837.33 


-1837.19 


tl 


389.102 


388.916 




-3.1742 


-5.3076 


ts 


11523.8 


11522.7 







40000 


J-St 





30000 


Xo 


0.411377 


0.410279 


Xi 


-0.832102 


-0.834832 


X2 


49.4875 


29.0669 


X3 


0.654962 


0.655322 


^! 





-0.5 







-3 


7 


0.3 


0.3 






2 






2 


Wo 


145.885 


146.048 


/+ 


1.000 


1.000 


fn 


1.044 


1.045 




1.055 


1.059 




1.050 


1.059 


Vw 


-2.00 


-2.06 


A 


320 


410 


V{v 


0.86 


084 


Ao 


28 


28 



Table 4. Rms (a) deviations between experimental data [T7] and HFB-17 or HFB17st 
predictions. The first line refers to all the 2149 measured masses M, the second to the 
masses Mnr of the subset of 185 neutron-rich nuclei with Sn < 5.0 MeV, the third to 
the 1988 measured neutron separation energies Sn and the fourth to 1868 measured 
beta-decay energies Qfs. The fifth line shows the comparison with the 782 measured 
charge radii [18] . Note that units for energy and radius are MeV and fm, respectively. 



HFB-17 HFB-17st 



cr (2 149 M) 



0.581 
0.729 
0.506 
0.583 
0.0300 



0.575 
0.738 
0.495 
0.585 
0.0302 



(T{Sn) 
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123456789 10 



r (fm) 

Figure 1. Top panel: neutron density of ^^Ca (solid line) and ^^Ca (dashed line) 
and spin-density in ^^Ca (squares). The densities are given in units of fm^''. Bottom 
panel: the same as top panel, but for ^°Ca and ^^Ca. 
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Figure 2. Difference between tfie nuclear mass of odd-A and odd-odd nuclei obtained 
with the BSkl7 force with and without the additional spin-density dependent terms. 
The calculation is made here assuming spherical symmetry. 
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Figure 3. Landau parameters for spin and spin-isospin channels as a function of 
density. Tlie left panel and middle panels show Gq and Gq, respectively, in symmetric 
nuclear matter and the right panel Gq in neutron matter. The SLySst and LNSst 
curves correspond to SLy5 and LNS standard Skyrme forces with the additional spin- 
density dependent terms, respectively [10]. The interactions BSkl7 and BSklTst are 
described in Sect. IV. 




Figure 4. Difference between the binding energies E/A{Ss,p) — E/A{Ss — 0,p) in 
symmetric matter for BSkl7 and BSkl7st Skyrme interaction as a function of the spin 
polarization Ss = (pt — Pi)/ P for the different densities, as indicated in the legend (in 
units of fm^'^). 




Figure 5. Difference between the binding energies E/A{ds,p) — E/A{Ss — 0,p) in 
symmetric matter for BSkl7st, LNSst and SLySst Skyrme interactions as a function 
of the spin polarizations and for p = 0.6 fm^'^. A different value from 0.0 to —3.0 for 
the parameter Xg is adopted for each line. 




Figure 6. Same as Fig. SI but for neutron matter. 



